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I. INTRODUCTION 

In his doctoral thesis R. P. Feynman, extending earlier 
work by Dirac [1], developed a novel formulation of quan- 
tum mechanics [2]. Unlike Schrodinger's formulation in 
terms of the solution of a partial differential equation 
Feynman based his description on the intuitive picture 
of paths connecting two points in space. Each path car- 
ries an amplitude and a phase that is given by its classical 
action. Incorporating the principle of superposition, the 
propagator, i.e. the probability amplitude to move from 
one point in space to another, is given by the sum (i.e. in- 
tegral) over all paths connecting these points. While this 
path integral formulation is equivalent to the standard 
Schrodinger theory, its implementation as an operational 
algorithm to solve quantum problems is complicated due 
the mathematical difficulties associated with the path in- 
tegration. One of its advantages is, however, the con- 
ceptual insight it can provide. Most importantly, mod- 
ern semiclassical theory invokes the convergence of the 
path manifold as contained in Feynman's path integral 
towards a discrete subset of classical paths of extremal 
action as ^ 0. 

The semiclassical approximation applicable at the bor- 
der between quantum and classical mechanics follows 
from the path integral formalism in the limit that the 
variation of the classical action is large compared to h for 
small path variations. Mesoscopic systems with linear di- 
mension D large compared to the de Broglie wavelength 
^dBi ^dB <^ represent prototypical cases for which 
semiclassical approximations are frequently invoked since 
in many cases ab-initio quantum calculations become un- 
feasible. Moreover, the description in terms of paths can 
provide detailed physical insights into spectral and trans- 
port properties. For example, dephasing and decohering 
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interactions are associated with a characteristic mean free 
pathlength ^mfp, thereby limiting phase coherent trans- 
port to short paths, £ < ^mfp- A well known example is 
the open quantum billiard in the ballistic regime which 
has been extensively studied in the last several years 
both experimentally and theoretically (see, e.g., Ref. 3-5 
and references therein) . The semiclassical approximation 
has contributed to the understanding of phase coherent 
transport effects such as conductance fluctuations and 
weal localization [6-26]. Most semiclassical approxima- 
tions to date have focused on either the spectral density 
p{E) [27] or transport coefficients determined for bal- 
listic transport by S matrix elements [10]. Semiclassi- 
cal calculations of the wavefunction ^(r,/c) itself which 
test the quantum to classical transition locally on the 
finest scale have remained a challenge. Few pioneering 
studies have been performed: scars in closed billiards 
could be reproduced by semiclassical calculations of the 
energy- averaged probability density (I ^(r, /c)p) [28]. For 
open chaotic billiards, statistical properties of the wave- 
functions like nodal point distributions have been found 
in good agreement with random wave models [29]. In 
the regime of high incident energies with a large number 
of open modes, quantum calculations [30] for the open 
chaotic stadium billiard have shown that the wavefunc- 
tion closely mirrors the path bundles [12] of short clas- 
sical scattering trajectories. A simple semiclassical ap- 
proximation to the wavefunction yields good quantiative 
agreement with the quantum wavefunction. 
We present in the following an accurate semiclassical de- 
termination of the fine-scale wavefunction ^(r, k) for an 
open ballistic billiard. We construct the wavefunction in 
terms of a sum over paths connecting the entrance lead 
with an arbitrary point r in the interior of the billiard 
closely following the Feynman path integral prescription. 
We employ the pseudopath semiclassical approximation 
(PSCA) [14-17] to include both classical and diffractive, 
i.e. non-classical paths into the path sum. We aim at 
a quantitative agreement with quantum wavefunctions 
for the low-energy regime with only few open modes in 
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the leads (quantum wires) and a semiclassical descrip- 
tion that pertains to the interior of the bihiard. We 
gauge the accuracy of the wavefunction by comparison 
with fuh quantum wavefunctions. For technical reasons, 
we focus on the rectangular billiard for which the enu- 
meration and summation of paths is still feasible since 
it is a prototypical example of a integrable system. We 
show that the convergence towards the quantum wave- 
function is controlled by the dwell time, or equivalent ly, 
by the mean pathlength of the scattering state at given 
energy. 

The outline of the paper is as follows: In Sec. II we briefly 
review the PSCA of the constant energy propagator and 
present its extension to calculations of wavefunctions of 
open quantum billiards. Concurrent calculations of the 
full quantum scattering state as well as its truncated form 
in which Fourier components associated with long path 
length are filtered out are discussed in Sec. III. In Sec. IV 
we present a quantitative comparison between semiclas- 
sical and quantum wavefunctions. We relate the conver- 
gence of the semiclassical wavefunction towards the exact 
quantum wavefunction to the dwell time, or, equivalently, 
to the Eisenbud-Wigner- Smith (EWS) delay time of the 
scattered wave inside the ballistic cavity. In the outlook 
(Sec. V) we briefly discuss future applications to leaky 
billiards and decoherence processes where long (coher- 
ent) paths or delay times are effectively suppressed and 
the PSCA offers a simple route to construct wavefunc- 
tions when short paths dominate the dynamics. 



coordinate system is the quantum point contact at the 
entrance lead denoted by ri = (x = 0, = 0) (see Fig. 1). 
The energy of the scattering state is given by E = 
with k = Y^/c^ -h {nir/dy. The potential in the interior 
of the billiard cavity vanishes, ^ = 0, and is infinitely 
high at the walls. Dephasing and decohering interactions 
with the environment are neglected. 
The classical dynamics of transport from the entrance to 
the exit lead (transmission) or back to the entrance (re- 
flection) is given by all classical paths p connecting the 
quantum point contacts (located at the center of the lead 
junctions) with each other. The action of a classical path 
p with path length Lp is given by Sp = kLp. The clas- 
sical phase space is structured in paths (or path bundles 
[12] when the finite size of the lead openings is taken into 
account) whose weight (area in phase space) is given by 
the deflection factor 



Dp{k) = 



dr'dr ' k drdk 
1 d^S 1 d^S 
k dr'dk ' /c^ dk^ 



(2) 



For the rectangular billiard the deflection factor is exactly 
Dp = l/{kLp). In Fig. 1 (b) we depict several classical 
paths connecting the entrance and exit point contact. 
Each path represents a path bundle and its color intensity 
is proportional to its weight, i.e. the deflection factor. For 
the rectangular billiard Dp is strictly positive and free of 
singularities. Consequently, the Maslov indev /ip entering 
semiclassics is given by 



II. SEMICLASSICAL THEORY FOR 
SCATTERING STATES 

A. The rectangular ballistic billiard 

We discuss in the following the PSCA for wavefunc- 
tions with the help of the specific example of a squared 
(rectangular) open billiard [Fig. 1 (a)]. The lead width 
d = 1/16 is small compared to the linear dimension D = 
1 in reduced units. Consequently, the short-wave limit 
XdB <^ D is reached for the internal dynamics of the cav- 
ity while the motion inside the leads (or quantum point 
contacts) is still in the quantum regime for low trans- 
verse mode numbers n with wavenumber kn = nir/d. All 
wavenumbers are given in the following in units of n/d. 
The PSCA is designed to describe the semiclassical dy- 
namics for billiards which are coupled to quantum wires. 
The asymptotic scattering boundary condition is defined 
by the incoming wave 

M/„,fe,(r) = ^Xn(y)e*'="^ (1) 

where (x, y) are the local coordinates along and 
perpendicular to the entrance lead and Xn{y) = 

sm[^{y — ^)] is the transverse wavefunction of 

mode n. ^^,^33(1*) is flux normalized. The center of our 



fip = 2Np (3) 

where Np is the number of reflections from the hard wall 
along path p. Note that one could alternatively incor- 
porate von Neumann boundary conditions by setting the 
Maslov index equal to zero. 



B. The PSCA for scattering states 

Starting point for the development of a semiclassical 
approximation to the wavefunction of the scattering state 
in the interior of the billiard is its expression in terms of 
the Green's function G(r', r, /c), 

nd/2 

vE^n(r', k) = -i / G(r^ x = 0, y)Xn{y)dy, (4) 

J-d/2 

for an incoming particle with wavenumber k = V2E in 
mode n. The scattering boundary condition (Eq. 1) at 
the lead entrance {x = O^y) acts as a source and the 
Green's function or constant-energy Feynman propagator 
connects the entrance point with the observation point 
via all (non-classical) paths. Eq. 4 contains an integral 
over all entrance points with transverse coordinate y. 
We note that the Fisher-Lee equation [31] expressing the 
S matrix elements in terms of G is a special case of Eq. 4. 




3 




FIG. 1. (Color online) (a) Geometry of the rectangular (square) billiard with side length D — 1 and equal lead widths d — 1/16. 
The placement of the leads is point symmetric with offset h — 0.25. The diffractive couplings from the leads (quantum point 
contacts) to the cavity, Cn{0,k), and for backscattering into the cavity v{0\0,k) are sketched. The quantum point contacts 
(red dots) placed in the center of the leads are located at the coordinates (xi = 0, 2/1 = 0) and (x2 — D^y2 — D — 2h). (b) The 
set of classical paths that connect the two quantum point contacts up to the maximal length L = 10. The color intensity of 
the paths is proportional to their deflection factor defined in Eq. 2. 



For example the transmission amplitude tmn is given by 

X y y x*M)G{x' = D,y' ,x = 0,y;k)xn{yW'^y, 

where the integral over y' extends over all coordinates at 
the exit lead weighted with the outgoing mode amplitude 
Xm{y') in the exit channel. The transmission (reflection) 
matrix is related to the S matrix as tmn = -S^n (^mn = 
S}^) and the conductance g of the quantum billiard is 
determined by the Landauer formula [32, 33]: 



9{k) 



2e^ 
h 



Op2 



N N 



(6) 



n=l m=l 



where N is the number of open modes in the leads (in 
our case the leads have equal width). 
The semiclassical approximation to Eq. 4 and Eq. 5 in- 
volves two steps: First, the evaluation of the integrals 
over the entrance/exit leads by either a stationary phase 
approximation (SPA) as done in the conventional semi- 
classical approximation or as a diffractive integral [12] 
as implemented in the PSCA. We employ for the lat- 
ter a combination of the geometric theory of diffraction 
[34] (GTD) and the uniform theory of diffraction [35, 36] 
(UTD), called GTD-UTD, which has been previously 
successfully applied for the scattering matrix of the cir- 
cular billiard [16]. The GTD-UTD takes into account 
the multiple scattering between the edges of a given lead 
(for details see [16]). Second, the quantum propagator 
is replaced within the PSCA by a semiclassical propa- 
gator that contains in addition to classical paths also 



pseudopaths, i.e. sequences of classical paths joined by 
diffractive (back) scattering at the lead openings. The 
reasoning underlying this augmented path manifold is 
that near the sharp edges of the leads or, more gener- 
ally, at the interface between the quantum point con- 
tacts and the cavity (Fig. 1), the semiclassical regime of 
sufficiently small A^^^ cannot be reached and, therefore, 
non-classical path contributions, determined to leading 
order in h by diffractive integrals, must be included from 
the outset. These pseudopaths are intuitively simple re- 
alizations of contributions to the Feynman path integral 
that are non-classical in origin, yet can be systematically 
included in the approximation. Accordingly, each pseu- 
dopath p of order 77 (p) consists of a sequence of 77 + 1 
classical path segments Pi joined by r] diffractive scat- 
terings at one of the lead openings (or point contacts) 
with diffractive amplitude v{9^\0^^~^ ^k)^ where 9^^~^ 
{0^'^) are the incoming (outgoing) scattering angle with 
which the trajectory pi_i (pi) approaches (leaves) the 
point contact. Explicit analytic expressions of the cou- 
pling coefficient v{9^%0^^~^ ^ k) are given in appendix A. 
The PSCA to the propagation along the pseudopath p 
between the starting point r and the end point reads 



G 



PSCA 
P 



(r',r,fc) = 



n 



-\k) 



(7) 



where the amplitude for each classical path segment pi 
connecting two quantum point contacts is given by the 
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standard SCA expression for the rectangular billiard 
27r 



(27ri)3/2 



Dp{k)\x 



(8) 



exp 



iSp{k) ^^Mp 



27T 



: exp [z/cLp — iNpir] 



The complete propagator G^^^^(r', r, /c) to order A is 
the sum over ah contributions from pseudopaths connect- 
ing and r with r]{p) < A: 

GfC^r',r,fc)= GfC^r',r,fc). (9) 

p:r]{p)<A 

Note that in Eq. 9 and r can be arbitrary points 
inside the bihiard. In the apphcation to the wavefunc- 
tion Eq. 4 we wih set r to be the entrance point contact 
ri and to be an arbitrary internal point inside the 
billiard while in Eq. 5 is one of the lead-billiard junc- 
tions when transmission or reflection is determined. The 
order A controls the degree to which diffractive contri- 
butions are included. Since the sum in Eq. 9 extends 
over infinitely many contributions its numerical evalua- 
tion requires in practice the limitation of path lengths 
by a cut-off length Lmax- For integrable billiards such as 
the rectangular billiard the number of trajectories below 
a maximum length M^^^(Lmax) increases quadratically, 
M^^^ (I/max) i^max- Howcvcr, the number of pseu- 
dopaths resulting from joining classical paths by a se- 
quence of diffractive couplings eventually proliferates ex- 
ponentially, MP^^'^(Lmax) cx exp(Lmax/^o). Therefore, 
sums over pseudopaths in numerical implementations can 
only by executed up to modest length Lmax- The con- 
vergence depends on the parameter pair (A, Lmax)- Note 
that Ga=o'^ (Eq. 9) is not equivalent to G^^^ (Eq. 8). 
In the SCA a path that hits the exit leaves the cavity. 
Within the PSCA the diffractive scattering at the exit 
gives rise to a plane wave and, to first order, a circular 
diffractive wave. In zeroth order the circular wave van- 
ishes but the plane wave leads to geometrically reflected 
paths which are not included in the SCA. 
The propagator within the PSCA, G^^^^, can be used to 
construct the semiclassical wavefunction via Eq. 4. How- 
ever, the integral over the entrance lead-billiard junc- 
tion weighted with the transverse mode wavefunction Xn 
would require the calculation of G^^^^{r\x = O^y) for 
all points —d/2 < y < d/2. For low mode numbers n we 
can approximate the integral by a diffraction approxima- 
tion that replaces the lead junction by a quantum point 
contact which acts as a point scatterer located at ri. 

We note that this description is appropriate for low- 
energy scattering with AdB ^ d while in the high energy 
limit AdB ^ d entire path bundles emanating from the 
finite-size lead opening rather than resolved paths conect- 
ing the point contacts should be included in the semiclas- 
sical approximation [30]. In the following we will focus 



on paths emitted from this point contact. The amplitude 
for a incoming particle in mode n to leave the point con- 
tact at the entrance lead with launching angle relative 
to the lead axis (see Fig. 1) and wavenumber k is de- 
noted by Cn{0, k) (see appendix A for its analytic form). 
Analogously cj^(6>, k) represents the amplitude for a tra- 
jectory incident on the junction under the angle to exit 
in mode m. 

Inserting Eq. 9 and Cn{0,k) into Eq. 4 yields the PSCA 
for the scattering wavefunction inside the billiard subject 
to the boundary condition (Eq. 1) of incident current in 
mode n 



^^^^{r\k)=-z^ (10) 
X Yl Gl^^^(v',v^,k)cn{ei,k), 

p:r]{p)<A 



where 0^ is the entrance angle of the trajectory p. Anal- 
ogously, the transmission amplitude follows as 



tmnik) — i\/kfYi kfi (H) 
X c™(^^,fc)GPSCA(r2,ri,A:)c„(0;,fc). 

p:r]{p)<A 

We will use in the following Eq. 11 as a complementary 
test for the accuracy of the PSCA to the wavefunction in 
Eq. 10. 

III. QUANTUM CALCULATIONS 

Before presenting results of the PSCA, we briefiy 
review the method employed for solving the underly- 
ing quantum problem. Our quantum calculations are 
based on the modular recursive Green's function method 
(MRGM) [37-39]. In the MRGM the two-dimensional 
(2D) Schodinger equation is solved numerically on a 
tight-binding grid, which leads to a non-quadratic (co- 
sine) dispersion relation, the main source of limita- 
tions within tight-binding approximations of the 2D 
Schrodinger equation. In the PSCA the dispersion re- 
lation has the correct form of E = k'^ /2. Thus, the com- 
parison between the PSCA and QM requires small grid 
spacings within the MRGM. We have used 30 grid points 
per half-wave length. A coarser discretization below this 
value leads to a visible shift in k of the QM transport re- 
sults with respect to the PSCA due to the non-quadratic 
dispersion relation. 

Even though the concept of paths and path lengths does 
not explicitly enter the quantum description, the scat- 
tering wavefunction ^(r, /c) can be Fourier analyzed in 
terms of its length component L, 



§(r,L) = j dk^{Y,k)e 



-ikL 



(12) 



As we will show the Fourier conjugate variable to k 
is closely related to the physical variable length of the 
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(semi) classical dynamics. For the numerical evaluation 
of Eq. 12, we perform a windowed Fourier transform in 

interval [/^min? ^max]- 

Eq. 12 is the generalization of the path length spec- 
troscopy [6, 12, 15, 40] of S matrix elements 

pkmax 

Smn{L)= / dkSmnik)e-"''^. (13) 

Smn {L) is the probability amplitude for a quantum path 
of length L to scatter from mode n to m. The path- 
length spectrum of open quantum systems decays for in- 
creasing length L. However, contributions from very long 
paths may become important near resonances of long- 
lived quasi-bound states. We will present examples of 
path length distributions entering ^ below. 
Assuming for the moment that the Fourier component 
^^(r,L) can be, indeed, identified with the semiclassical 
path length, it is now instructive to construct truncated 
quantum wavefunctions that retain only Fourier compo- 
nents with L less than the maximum path length L^ax 
included in the PSCA. To this end, we truncated the 
inverse Fourier transform at L = L^ax 

^^{r,k)= dL §(r,L)e^^^. (14) 

Jo 

Analogous truncation of Fourier spectra are performed 
for S matrix elements 

Sl„{k) = / dL 5™„(L)e*^ (15) 

Normalizations have been omitted for simplicity. The 
comparison between and ^psca g^^^Q^g directly 
and quantitatively compare the semiclassical wavefunc- 
tion with the quantum wavefunction that contains all 
Feynman paths up to the same length Lmax- Conversely, 
comparison between and the full quantum wave- 
function ^ allows to assess the infiuence of long paths 
L > Lmax and, therefore, the truncation error involved 
in semiclassical path sums (see Sec. IV B). 
The numerical evaluation of the Fourier transform for 
finite discretized intervals in k gives rise to a maximal 
resolvable length AL = 27r/dk where Sk is the grid spac- 
ing in the k domain. In order for contributions with L > 
Lmax not to enter the amplitude at the maximum re- 
solvable length AL, ^(r, AL), must already be strongly 
suppressed. Otherwise, the Fourier spectrum is back- 
folded such that contributions with L > AL ^ L^ax 
appear near the origin (L ^ 0) and cannot be cleanly 
cut off by Eq. 14. To avoid such back-folding we choose 
a large AL. However, we will see in Sec. IV B that a 
complete truncation cannot be established in the vicin- 
ity of sharp resonances. This truncation is known as the 
sine filter in signal processing. It leads to a smoothing 
of sharp peaks on the k scale and causes a violation of 
unitarity since high Fourier components corresponding to 
contributions from long paths are missing. We note that 



the energy average of the local density of states in the 
semiclassical description of scars [28] corresponds to an 
implicit truncation scheme. 

The relative importance of long paths can be quantified 
by the expectation value of the Eisenbud-Wigner- Smith 
(EWS) time delay operator for incoming mode n 

{Q) = {-iS^^S) 

m 

For billiards with zero potential in the interior, the time 
delay tews can be directly converted into a path length 

^EWS kTEWS' (17) 

By comparing ^ews(^) at at given wavenumber (or en- 
ergy) with Lmax we can provide an independent estimate 
for the expected proximity of to ^ and, in turn, the 
convergence of ^p^ca ^^^i scattering state. 

An alternative measure for the time the particle spends 
inside the cavity is the dwell time 

rD= J \^{r,k)\^dxdy. (18) 

cavity 

The difference between tews and td is referred to as 
the interference delay that the wave packet experiences 
before entering the cavity due to interference with parts 
of itself that have already been refiected [41]. This 
self-interference delay becomes important when the de 
Broglie wavelength of the particle is comparable to the 
linear dimension of the billiard. In the semiclassical 
regime, however, we find td = tews to a very good 
degree of approximation. 



IV. NUMERICAL RESULTS 

A. Scattering matrix 

To set the stage, we first present typical results for S 
matrix elements (Eq. 5 and Eq. 11). The total transmis- 
sion T{k) (Eq. 6) as a function of the incident energy 
or, equivalently, k (Fig. 2) displays excellent agreement 
between the PSCA and the scattering matrix trun- 
cated at the same path length L^ax = 17.5 as the PSCA. 
A few exceptions are worth mentioning. They appear, 
e.g, in the vicinity of the channel opening k < 27T/d. 
Just above the channel opening the emission angle is 
close to 7r/2. For grazing incidence the diffraction ampli- 
tude Cn{0^ k) and scattering amplitude v{9' ^ 0^ k) in GTD- 
UTD is less accurate, most likely causing this discrep- 
ancy. While overall the agreement between the PSCA 
and the truncated quantum S matrix is remarkable, pro- 
nounced differences appear to the full quantum S ma- 
trix as highlighted in Fig. 2 for the magnified interval 
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k in Tv/d 

FIG. 2. (Color online) Comparison of the total transmission T(k) as a function of k within the PSCA (green line), truncated 
QM (black dashed line) and full QM (in the inset, upper panel red line). The PSCA is calculated with a maximum path 
length of Lmax = 17.5 and maximum order of diffract ive scattering A = 6. The inset shows a magnification of the region 
k e [1.75,1.85]. 



1.75 < k < 1.85. Pronounced sharp structures are miss- 
ing in the PSCA indicating that the latter are due to 
long paths, i.e. represent long-lived resonances. In prin- 
ciple, the PSCA could account for those if Lmax could be 
extended. In practice, however, the exponential prolif- 
eration of contributing pseudopaths prevents to perform 
complete path sums. The comparison suggests that for 
those k values where sharp resonances appear, the PSCA 
wavefunction will be significantly different from the quan- 
tum scattering state ^ while in spectral regions where 
T{k) is smooth, the PSCA should become accurate. 



B. Wavef unctions 

We present in the following three prototypical cases of 
scattering wavefunctions: one in the smooth off-resonant 
part of the spectrum for incoming mode n = 1 (Fig. 3), 
one near a sharp resonance (Fig. 4) and one for incom- 
ing mode n = 2 (see Fig. 10 below). Obviously, in the 
non-resonant case we find near perfect agreement be- 
tween ^, and ^PSCA [^i the resonant case 
^PSCA gtrongly differs from ^ but is in close agreement 
with as anticipated. In resonant scattering a long- 
lived quasi-bound state is excited such that contributions 
from very long path lengths become essential while off- 
resonant wavefunctions typically show signatures associ- 
ated with the dynamics of short paths. 
The computational effort to calculate ^psca much 
larger than for the S matrix elements since the path sum 
has to be performed for each grid point of a square lat- 
tice with small spacings Sx. In order to generate high- 
resolution images of ^p^ca ^[^i^ modest computational 
effort we employ the following trick: we expand ^p^ca 
in terms of the analytically known wavefunctions of the 



closed billiard 



{r\mn) — sin (Kmx) sin {KnV), 



(19) 



where Kn = With the help of the projection ampli- 
tude 



.PSCA 



(20) 



evaluated on the grid, we can evaluate ^^^ca 
terior as 



f PSCA(r, k) = -J2 a^Jn^ik) sin (K^x) sin (K^y). 

ran 

(21) 



Due to energy conservation the expansion coefficients 
o^n^{^) non-zero only near the circle K'^^K'^ = k'^ 
(see, e.g., the inset of Fig. 5). In the limit of a bound 
state of a closed billiard the amplitudes take the form 
of a delta function amn{k) = S{k — \JK^ + ^m)- We 
can thus use a coarse grid 5x <T:/k such that the max- 
imal resolved wavenumber Ai^ = / 5x is slightly larger 
then the wavenumber k of the scattering state. amn{k) 
calculated from this coarse grid contains the complete 
information on the wavefunction inside the cavity. We 
note that the eigenfunctions of the closed billiard do not 
form a complete basis for the open billiard because the 
lead openings are replaced by hard walls. Therefore, the 
Fourier expansion (Eq. 20) and its inverse are not strictly 
unitary. 

Since the expansion coefficients amn are non-zero only 
near the circle k'^ = K^-\-K^^ we can construct a spectral 
angular amplitude of the scattering state with wavenum- 
ber k as 

Q{(j)^k) = ^^^5{n^tdin{(j))m)amn{k). (22) 
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FIG. 3. (Color online) Comparison of x^f^^^^^ the trun- 
cated wavefunction and the full wavefunction ^ for a 
non-resonant scattering state in incident mode n = 1 with 
k — 1.7835 and ^ews = 10.5. (a) Full quantum wavefunction 
^, (b) truncated quantum wavefunction with Fourier compo- 
nents L < Lmax = 17.5, (c) PSCA with Lmax = 17.5 with 
maximum diffractive order A = 6. 



FIG. 4. (Color online) As Fig. 3 but for a resonant scattering 
state in incident mode n = 1 with k = 1.8295 and fews = 47. 
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The expression (5(n, tan(^)m) is unity for tan(^) = n/m 
and zero elsewhere. Eq. 22 can be apphed to both the 
exact quantum state and to its semiclassical approxima- 
tion. Fig. 5 shows the angular spectra of the non-resonant 
scattering state [Fig. 3 (a)]. As expected from the visual 
inspection of the wavefunction, the angular distribution 
shows many peaks corresponding to a large number of 
excited modes. Nevertheless ^{(j)^ k) is not uniformly dis- 
tributed in angle (j) because excitation of modes with large 
(j) would lead, in general, to long lifetimes caused by their 
weak cavity-lead coupling [see Fig. 9(a)]. For states with 
short dwell times such as Fig. 3 large (j) contributions are 
suppressed. 

To quantify the agreement between the quantum and 
semiclassical wavefunctions we use the distance metric in 
Hilbert space 



d'(*l,*2) 



1 



dxdy |«'2(r,fc) - «'i(r,fc)p. (23) 



cavity 



The convergence of the semiclassical wavefunction ^psca 
as a function of A and Lrnay. 

can be conveniently studied 
by employing ^^^^^). We observe that at fixed 

Lmax and increasing A, ^p^ca converges monotonically 
to which is shown in Fig. 6 (a) for the scattering state 
at /c = 1.7835. The semiclassical wavefunction ^psca 
has converged to within dP{^^ , ^psca^ ^ 0.0015 for 
A = 6 at L^ax = 20. 



0.025 




FIG. 5. (Color online) The absolute square of the spec- 
tral angular amplitude ^^(0, k) (Eq. 22) at fixed wavenumber 
k = 1.7835 for ^ [Fig. 3 (a)]. The inset shows the expansion 
coefficients anm{k) and the definition of the angle 0. 



At fixed A and increasing Lmax the semiclassical wave- 
function ^PSCA gi^g^^^i^g diverge from with the on- 
set of divergence shifted to larger Lmax as A increases 
[Fig. 6(a)]. This is a consequence of the increasing lack 
of pseudopaths required for the complete path sum at 
length Lmax- For large L^ax an exponentially increasing 
number of pseudopaths of the same length become acces- 
sible that consist of segments of shorter classical paths 
joined by an increasing number of diffr active scatterings 
(see Eq. 7). Some of these high-order pseudopaths are 
missing for the proper interference with the paths already 
included in the PSCA of low order A. Therefore, conver- 
gence to the quantum wavefunction requires a correlated 
limit of both Lmax ^ oo and A ^ oo. It should be 
noted that with increasing order A, the PSCA becomes 
more and more sensitive to the accuracy of v{9' ^ 9^ k) (the 
diffractive amplitude for the internal scattering at the 
cavity-lead junction) since any error in v{6\6,k) is ex- 
ponentiated to the power A. 

While the dependence of ^psca^ order A 

of the PSCA allows to estimate the significance of diffrac- 
tive contributions to the path sum up to a given length 
^max, the distance (i^(^,^^) measures the total contri- 
bution of paths with length beyond L^ax [Fig. 6(b)]. The 
convergence of to the exact scattering state as a func- 
tion of Lmax is controlled by the mean path length ^ews 
of the scattering state. In general, for L^ax ^ ^ews the 
truncated wavefunction and thereby ^psca g^£_ 
ficiently high A converge to the full scattering state. 
The PSCA fails to reproduce the quantum scattering 
state near a sharp resonance (Fig. 4) because long paths 
well beyond Lmax contribute to this long-lived quasi- 
bound state. It is therefore instructive to directly com- 
pare the distance functions (i^(^, ^p^^^) with the EWS 
length ^Ews as a function of k (Fig. 7). Indeed, the 
distance of the wavefunction in Hilbert space strongly 
correlates with the mean path length ^ews and allows 
to predict the accuracy of ^p^ca ^ given L^ax when 
^Ews is known. The two prototypical cases shown above 
are marked in Fig. 7. It is also instructive to measure 
the distance between ^psca (^Yig. 7). It is uni- 

formly small (< 0.05) over the entire k interval and slowly 
varying. It is worth noting that the residual fiuctuations 
in (i^(^^, ^^^^^) are not due to shortcomings of the 
PSCA but rather due to the truncation process in (see 
Eq. 14). Near resonances ^(r, L) contains non-negligible 
contributions up to very long paths. Consequently, when 
^(r, L) is still non- negligible for L beyond the Fourier 
resolution limit L > AL = 27r/^/c, back folding causes 
truncation errors. The latter are responsible for the in- 
crease of (i^(^^, ^P^^^) in the region of high density of 
sharp resonances (e.g. near k = 1.765 and k = 1.822). 
The contribution from long paths and the possible con- 
tamination by back folding can be conveniently moni- 
tored by the projection amplitudes amn{k) onto bound 
states of the closed billiard (Eq. 20) or onto quasi-bound 
states of the open billiard. We present |a^n(/c)p and its 
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FIG. 6. (Color online) (a) The distance as a function of Lmax depicted for varying orders A = 1 to A = 6 and 

wavenumber k = 1.7835 (same as in Fig. 3). (b) The normalized distance ^2(^'o) of ^ to the exact quantum wavefunction ^ 

d2(^'o) with ^ being either ^^^^^^ (green line) or (black line) as a function of Lmax for the same wavenumber as in (a). 
The mean path length within the full QM (Eq. 17) is ^ews = 10.5 and is marked by a vertical line. The PSCA is depicted 
for order A = 6. (black line) is obtained via a discrete Fourier transform of Eq. 14 in the interval k G [1.725, 1.875] with 
6k = 0.00025. 



path length spectrum |arnn(^)P, where 

dmn{L) = I dk a^n{k)e-'^'^ (24) 

for two extreme cases in Fig. 8, the broad resonances 
(m = 29, n = 1) and the sharp resonances (m = 15, 
n = 25). For the broad resonance, PSCA and the trun- 
cated quantum state agree perfectly since the path length 
spectrum of the amplitude has already decayed by over 
four orders of magnitude near L = AL = 500. For the 
sharp resonance with ^ews = 1460, the path length spec- 
tral intensity at L ^ 500 is still ~ 10~^ of its value for 
small L and, consequently, truncation leads to large dis- 
crepancies between ^psca Naturally both the 
PSCA and the quantum state truncated at L^ax = 17.50 
fail to describe the exact quantum scattering state for 
the sharp resonance. 

The proximity of the scattering wavefunction to a bound 
state of the closed billiard can also be monitored by the 
angular spectrum (Eq. 22). The sharp resonance de- 
picted in Fig. 4 is very close to (m = 4, n = 29). Con- 
sequently, a single peak strongly determines the angu- 
lar spectrum at tan^ = n/m overshadowing all other 
contributions [see Fig. 9 (a), note the different scale in 



the inset]. For the broad resonance in Fig. 3 the angu- 
lar spectrum contains many components of comparable 
magnitude [see Fig. 9 (b)]. In this case, 
ft closely agree with each other. 

While the previous numerical examples feature low 
wavenumbers k < 2 with only one open channel we 
have checked on the convergence of ^psca higher 
wavenumbers up to A: = 5 as well. We find, in general, 
the same very good agreement between the PSCA and 
the full quantum mechanics in the regime L^ax > ^Ews 
and between the PSCA and the truncated quantum me- 
chanics in the regime L^ax < ^ews- 
An example for the wavefunction convergence for the sec- 
ond mode with k = 2.25375 and ^ews = 9.5 is shown 
in Fig. 10. One remarkable feature of this wavefunction 
different from the previous cases is the intensity enhance- 
ment along a classical periodic orbit (see Fig. 10). Both 
^PSCA g^^^ g]^Q^ clear traces of this orbit. For the 
full quantum state ^ the occurrence of this structure is 
caused by the excitation with nearly equal amplitude of 
two almost degenerate eigenstates of the closed system 
\mi^ni) with mi = 26, ni = 25 and m2 = 30, n2 = 20. 
Superposition of these two eigenstates lead to an enve- 
lope that follows the track of the periodic orbit similar 
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k = 1.82957r/d 




k in Tv/d 



between the quantum wavefunctions ^ and the PSCA ^^^^^ [red (dark gray) 



FIG. 7. (Color onhne) Distance d'^{^,^^' 

hne], as well as the distance (i^(^^, ^^^^^) between the truncated QM and the PSCA [green (light gray)line]. The PSCA is 
calculated in order A = 6 and with Lmax = 17.5. The black solid line is the mean path length ^ews (see Eq. 17) and the 
horizontal dashed line marks a mean path length of ^ = i^max = 17.5. Note the logarithmic scale of the figure. The two 
vertical arrows mark, respectively, the wavenumber of a non-resonant scattering state (k — 1.7835, same as in Fig. 3) and the 
wavenumber of a resonant scattering state {k — 1.8295, Fig. 4). The sharp peak oX k — 1.82925 corresponds to the resonance 
m = 29, n = 4 and has a mean path length of ^ews = 145 almost degenerate to the resonance state m = 4, n = 29 at A; = 1.8295. 
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FIG. 8. (Color online) The absolute magnitude square of the amplitudes \amn{k)\^ for the two extreme cases of (a) a broad 
resonance (m = 1, n = 29) and (b) a sharp resonance (m = 25, n = 15). Note the different ranges of the ordinate in (a) and (b) 
for the quantum result. The inset shows the length spectrum \amn{L)\^ of the amplitude of ^ in a logarithmic scale. The length 
spectrum of the sharp resonance [inset of (b)] is not sufficiently decreased at the boundary of numerical resolution AL = 500 
such that some paths with L > Lmax are not filtered out of This causes the difference between (black dashed line) 
and ^P^*^^ (green line) in (b). The scattering state which excites the sharp resonance corresponds to the mean path length of 
£ = 1460 at /c = 1.822 (see Fig. 7). The PSCA is calculated with Lmax = 17.5 and order A = 6. 



to the formation of a beat in acoustics [42] . It is now in- 
structive to analyse their appearance within the PSCA. 
We first note that the angular spectrum of ^ (analo- 
gous to Fig. 9) features two dominant peaks with equal 
height near the angles = tan~-^(ni/m^) {i = 1,2) de- 
termined by the quantum numbers of the two resonant 
states. Moreover the angular distribution for diffractive 
coupling from the point contact into the billiard in sec- 
ond mode C2(^, k) (Eq. 10) features a high amplitude near 



these angles, 0^ = This differs from the coupling in 
first mode for which large angles would be suppressed. 
The build-up of the beat pattern results now from the 
constructive interference between the rays emitted from 
the contact into the billiard with entrance angle close to 
and its replica propagating into the same direction af- 
ter 4 consecutive refiections at the billiard wall. The path 
length difference between two such trajectory inside the 
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FIG. 9. (Color online) Comparison of ^^(0, k) within PSCA, truncated QM and full QM for a non-resonant (k = 1.7835, Fig. 3) 
and a resonant {k = 1.8295, Fig. 4) scattering state. The wavenumbers correspond to those labelled in Fig. 7. Lmax = 17.5, 
A = 6. 




(b) PSCA (c) truncated QM 






FIG. 10. (Color online) As Fig. 3 but for an off-resonant scattering state in incident mode n = 2 with k = 2.25375 and 
^Ews = 9.5. Comparison of ^^^^^^^ the truncated wavefunction and ^. The white dashed line in (a) shows the density 
enhancement along the classical periodic orbit. 



billiard is given by 



2D 



rrii + Hi 



(25) 



The necessary condition for constructive interference 

AS' = kAL = 27rj (j G N) (26) 



is now approximately satisfied for both (f)i {kAL = 
IOO.OItt) and 02 {kAL = 101.977r) at = 2.25375. In 
addition the condition for constructive interference be- 
tween the rays with entrance angle 0^ = 0^ and 0^ = 
leads to sm{hky) = 0, which gives a primitive estimate for 
the strength of the coupling between eigenstate and in- 
coming second mode (for definition of h see Fig. 1). Both 
requirements are approximately met for (m^,n^)i^i^2 re- 
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suiting in the simultaneous excitation of the correspond- 
ing eigenstates and, in turn, in the beat pattern and in 
the density enhancement near the periodic orbit. The 
strong coupling between the incident mode n and the 
two eigenstates is facilitated by the spacial proximity, 
to within a de-Broglie wavelength, of one of its impact 
points at the wall to the location of the entrance lead. 
The excitation amplitude of the eigenstate also depends 
on the position of the exit lead. Taking into account 
higher-order diffractive corrections incorporates the cou- 
pling to the exit lead into the semiclassical description. 
Note that the present mechanism for the formation of this 
periodic orbit is different from the appearance of scars 
in wavefunctions of open chaotic billiards for low k [43] 
which is based on constructive interference for consecu- 
tive retracing of a single isolated periodic orbit. It is also 
different from the appearance of path bundles for large 
k where bundles of short classical scattering trajectories 
emanating from the entrance opening whose width is de- 
termined by the lead width appear in scattering states 
[30, 44]. 

The comparison between the truncated and the full QM 
wavefunction (Fig. 10) shows that there are non- vanishing 
contributions from paths with length L > L^ax indicat- 
ing the presence of resonances in the immediate vicin- 
ity. Investigation of as function of the trunca- 
tion length Lmax shows rapid decay for small L^ax and 
a slow, plateau-like, decay for L^ax > H due to the 
residual influence of nearby resonances which causes a 
small ^Ews but still non- vanishing distance at 

^max = 17.5 > ^EWS- 



V. SUMMARY AND OUTLOOK 

We have presented the construction of the semiclassi- 
cal constant energy propagator and of scattering states 
employing the pseudopath semiclassical approximation 
(PSCA). The convergence of the PSCA to the quantum 
limit is controlled by the maximum path length L^ax and 
the maximum order A of non-classical diffractive scatter- 
ings included. For the open rectangular billiard we find 
unprecedented quantitative agreement between ^psca 
and the full quantum scattering state ^ when the mean 
path length ^ews = ^''"ews determined by the Eisenbud- 
Wigner-Smith (EWS) time delay tews is covered by the 
PSCA, ^Ews ^ ^max- Thus, the pseudopaths resulting 
from sequences of classical paths joined by non-classical 
diffractive scatterings at the lead openings provide the 
necessary complements to the classical paths for comple- 
tion of the Feynman path sum for quantum propagation. 
To our knowledge, this is the first protocol for construct- 
ing semiclassical scattering wavefunctions whose conver- 
gence to its quantum counterpart can be quantitatively 
controlled. Even when the mean path length of the exact 
scattering state ^ews exceeds the maximum path length 
included in the numerical implementation of the PSCA 
as it happens for energies (or wavenumbers k) near long- 



lived resonances, we find near-perfect agreement between 
^PSCA g^^^ corresponding quantum wavefunction 
with the path length spectrum truncated at the same 
^max as the PSCA. We conclude by pointing to a few 
future applications and extensions: due to the math- 
ematical equivalence between the Schrodinger and the 
Helmholtz equations the experimental measurement of 
such wavefunctions can be conducted in open microwave 
billiards using movable antennas [45] or in micro cavity 
lasers [42]. The direct measurement of wavefunctions in 
quantum dots is still a major challenge but some progress 
has been made using scanning tunneling microscopy in 
graphene quantum dots [46]. The truncation of long 
paths beyond L^ax, introduced here to control the sum 
over an exponentially proliferating set of pseudopaths, 
does have, in fact, experimental analogues and applica- 
tions: finite energy resolution of the detection and/or 
excitation processes leading to smearing out of the sharp 
resonances is equivalent to suppressing long paths in the 
expectation value ([^(r)^)^;. Moreover, the finite phase 
coherence length present in decohering systems pre- 
cludes the appearance of long-lived resonances and causes 
long paths to contribute only incoherently [7, 47, 48]. The 
latter can easily be incorporated within the PSCA by an 
exponential damping e~^/^^ of long path contributions 
[6]. Damping of long paths naturally occurs in billiards 
with leaky boundary conditions where tunneling through 
the billiard walls is present. In the experiment such sys- 
tems can be realized as microwave billiards with dielec- 
tric boundaries. The PSCA can be used to describe such 
systems by incorporating an additional reflection ampli- 
tude for each bounce off the billiard walls accounting for 
the tunneling probability. For such billiards, violation of 
flux conservation (unitarity) as well as a small mean path 
length, both naturally incorporated within the PSCA, 
are key features of the scattering system. The appear- 
ance of total reflection for small incident angles causes 
the reflection amplitude to become a pure phase factor, 
i.e. its modulus is unity. In the ray picture the asso- 
ciated phase shift can be interpreted by a spatial shift 
known as the Goos-Hanchen shift [49, 50]. The inclu- 
sion of this effect can lead to a considerable change in 
the dynamics of billiards with penetrable walls [51-54]. 
Due to the possibility of controlling individual path con- 
tributions, the PSCA has the potential to develop into 
an accurate method for calculating scattering states in 
open billiards where both diffraction at lead edges and 
dielectric boundary conditions are present. 
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Appendix A: Diffraction amplitudes 

We reproduce here the diffraction amphtudes v{6\ 6^ k) 
and CmiO^k) within the GTD-UTD. The derivation can 
be found in Ref. 16. 

The diffraction ampUtude v^™{0\6,k) for backscatter- 
ing into the cavity within the GTD [34] is given by 



.,GTD//)/ 



1 



{e',e,k,d) = -DL{e',e)e 



ifef (sin6>'+sin6') 



(Al) 

with the scattering coefficients at the left and right wedge 

DL{e',e) = D{^/2-e\^/2-e) 

DR{e',0)=D{TT/2 + 0',TT/2 + e) (A2) 

and 



D( 



sin7r/7V 



1 



cos 



N 



COS 



N 



COS 



N 



COS 



N 



(A3) 



with N = 3/2 the exterior angle (in units of tt) of a 
perpendicular wedge. The angles 6 and 6^ are depicted in 



Fig. 1. We use the diffraction coefficient within the UTD 
[35] to take into account multiple scatterings between the 
two wedges of the cavity- lead junction: 



D 



UTD / // 



(t),r',r,k) 



N 



(J cot I ^ 1 X 



cr, 77 = ±1 

F I k 



2N 



(A4) 



where a±{j3) = 2cos^ ^t^Nu^ ^± ^j^^ integer 

which most closely satisfies 27i Nn^ — P = i^r. The func- 
tion F is defined as a generalized Fresnel integral: 



F{x) = -2^^fxe-''' f 
J X/ 

With the notation 



dre' 



UL{e', e, r, k) = L>u™(7r/2 - 9', -k/2 - 0, r' 
Ur{0', 0, r, k) = Z)"™(7r/2 + 0' , 7r/2 + 0, r' 



(A5) 



oo, r, k) 
oo, r, k) 
(A6) 



we obtain for the diffraction amplitude viO'^O^k) within 
the GTD-UTD: 



. Jmax 

- -^/2, j^, k)gj{k)e'^~^DR{^^/2, 0) + /7^(^^ ^^/2,jd, k)gj{k)e^^^~ Dl{-^ /2, 6) + 



odd: j = l 

Jmax 



-| jmax 

T Jd, %,(fc)e**^^Z)fl(+7r/2, 0) + Ul{0', -7r/2, jd, k)g,{k)e'^" Dl{-^ /2, 0) 



even: j = l 



(A7) 



where 



and 



,i{kjd+{j -1)7,) 



A:^(±sinl9'± sin l9). 



(A8) 



(A9) 



The sum goes over multiples of scatterings between the 
wedges and is cut at jmax where convergence is reached. 

We use jmax = 5. 

The diffraction amplitude for coupling of the quantum 
lead to a cavity within the GTD is given by 



{0,k,d) 



xm-n 



2dk'y 



-DLi0,0m)e''^e 



-ik^ sin 9 



1 



gi/c ^ sin 9 



(AlO) 



As before we use the notation 



~2 
37r 

Y 



(All) 



for the left and right wedge. The angle dra is deter- 
mined by the open lead mode and is given by dm — 
arcsin(m7r/(i/c). Using the UTD Eq. A6 for the out- 
coupling into the cavity we obtain: 
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OTD-UTD 



{0,k,d)=c^^''{0,k,d)-i 



e 2 



J2 UR{0,+Tr/2,jd,k)gj{k)e"t'^^DL{-7r/2,0^) - -7r/2, jd, %,(fc)e^*""Dfi(+7r/2, 0^) 

Lodd:j=l 



£ Unie, +n/2Jd, k)gj{k)e'^^~DR{+n/2, 0^) + Ul{0, -7r/2, jd, k)gj{k)e'^~* Dl{-^ /2, 9^) 

even:j = l 



(A12) 



where 



and gj{k) is given in Eq. (A8). 



= ±k- smO ± , 

2 2 ' 



(A13) 
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